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1. Introduction 

Numerous investigations of nonequilibrium physics have been carried out recently, 
and some general results have already been found. For example, the fluctuation 
theorem jTHTQ], the Jarzynski equality [TT}[T2], an< ^ the Hatano-Sasa relation [13] are 
examples of recent developments in nonequilibrium physics. These concepts are related 
to each other [I1HT6] . and then one could imagine that they may include interesting 
mathematical structures. Moreover, if there are mathematically interesting structures, 
it could be expected that these mathematical structures are also physically-meaningful. 
Hence, it would be valuable to seek mathematical structures behind nonequilibrium 
states. 

We here focus on stochastic processes described by master equations, and try to 
find mathematical structures behind them. In general, a master equation describes 
the time-development of the probability of states. One of the important quantities 
for nonequilibrium states is the 'current', which is related to the number of specific 
transitions occurring during a given time. For example, in a particle hopping model 
on a one dimensional lattice, the numbers of particles moving in and out of a specific 
site gives information about the current. In addition, in order to evaluate all statistics, 
including fluctuations, of the number of specific transitions, the framework of counting 
statistics is available [T71 - [2Tj . Actually, the counting statistics for stochastic processes 
enables one to obtain information about nonequilibrium currents in a simple two-state 
model under cyclic perturbations, for both adiabatic regimes [22TI25] and nonadiabatic 
regimes [26]|2T] ; here, 'adiabatic' means that the perturbations are performed very slowly. 
In general, a cumulant generating function for the counting statistics can be interpreted 
as a 'phase', as discussed in [22]. The phase is divided into two parts, i.e., the dynamical 
phase and a remaining one. In previous studies [22H27] , it has been shown that the 
concept of geometric phase, such as the Berry phase and the Aharonov-Anandan phase, 
is useful to calculate the current statistics due to cyclic perturbations. 

While a noncyclic geometric phase has already been studied in quantum mechanics 
[28H31] . the geometric phase discussion in the counting statistics has been extended to 
noncyclic cases recently [32j|33]. However, the discussions in [32j[33] are restricted to 
adiabatic cases, which correspond to quasi-static evolutions. 

The aim of the present paper is to develop a geometric discussion for noncyclic 
and nonadiabatic cases. In order to investigate transitions between two nonequilibrium 
states, one has to consider nonadiabatic and noncyclic evolutions. In the present 
paper, we will show that the remaining phase has a geometric interpretation even 
in the nonadiabatic and noncyclic cases, based on the discussions given by Samuel 
and Bhandari for quantum systems [28]. Hence, the present paper will give a natural 
extension of the adiabatic case given in [32}[33] to the nonadiabatic case. 

The present paper is constructed as follows. In section 2 we briefly explain the 
general scheme for counting statistics, and rewrite the formulation in terms of bra and 
ket state vectors. Section 3 is a brief explanation for cyclic evolution cases. Section 4 is 
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the main part of the present paper; we confirm that the remaining phase in noncyclic 
cases is adequately interpreted as a geometric phase. In section 5, we give concluding 
remarks and some discussions about physical meanings of the noncyclic geometric phase. 



2. Time evolution equation for counting statistics 

2.1. Basics for counting statistics 

We here briefly review counting statistics for stochastic processes. The formulation is 
based on a usual generating function approach, which has been widely used in various 
contexts (e.g., for the charge transport or full counting statistics in quantum systems, 
see [53].) Details concerning the counting statistics can be found in the Appendix or, 
for example, in [19] . 

Let us start from a master equation, which describes the time evolution of the 
probability distribution. The master equation is given by 

T7Pn(t) =Y, K nm(t) P m(t), (1) 
at m 

where p n {t) is the probability of finding the system in state n, and K nm (t) is the rate 
constant for transition m — > n. In general, the transition matrix {K nm (t)} is time- 
dependent. 

The aim of the counting statistics is to count the numbers of target transitions. For 
simplicity, we here consider a case in which the number of target transitions i^ — > Ja is 
counted. It is easy to extend the following discussions to multiple transition cases, and 
in these cases it is possible to derive a joint probability distribution for several target 
transitions. 

In order to extract statistics for the target transition i^ — » ja, a generating function 
is defined as 

oo 

F(x,t) = J2e k ^P(k\t), (2) 

k=0 

where P(k\i) is the probability of k target transitions za — > Ja during the time t. 
Derivatives of F(x, t) with respect to x & ve various information about the statistics 
of the target transition. As shown in Appendix, the generating function F(x, t) is 
constructed from a linear combination of functions {f n (x,t)} : 

F(x,t) = J2fn(x,t). (3) 

n 

fn(Xft) is a restricted generating function under the condition that the system is in 
state n at time t, and it satisfies the following time-evolution equation: 

~nfn{X,t) =Yl K nm{X,t)f m {x,t). (4) 

at m 

The matrix {n nm {Xi t)} is a modified transition matrix obtained by introducing counting 
parameters x m to the transition matrix {K nm (t)}; the off-diagonal element Kj AiA (t), 
which is the rate constant for the target transition, is replaced by Kj AiA (t)e x [35]. Note 
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that the modified transition matrix becomes the same as the transition matrix if we set 
X = 0, i.e., {n n m(x — 0,£)} = {K nm (t)}. In addition, (TjJ should be solved using the 
initial condition f n (Xify — Pn(0) (see Appendix). 



2.2. Expressions in terms of bra and ket state vectors 

The time-evolution equations for the restricted generating functions (j4|) can be rewritten 
in terms of a 'ket' state vector: 

^mx,t)) = H{ X ,t)\i(>{x,t)), (5) 

where H(x,t) = {ft nm (x, t)} is the modified transition matrix. The ket state vector 
\ip(x,t)) is defined as \tp(x,t)} = (fi(x, t), h{x,t), ■ ■ -) T ; where the superscript (-) T 
denotes the transposed matrix. 

Since all solutions to (jHJ) form a linear vector space, there is a dual space. In fact, 
we can construct the adjoint equation of ^ as follows: 

^(Hx,t)\ = -(iP(x,t)\H(x,t). (6) 

Note that (i/j(x, t)\ is not the complex conjugate of \t/)(x,t)), because of the non- 
Hermitian property of H(x,t). In addition, we can show that the norm, i.e, the inner 
product between (i/;(x,t)\ and \i/;(x,t)), is preserved: 

f t m X ,t)mx,t))=0. (7) 

This conserved norm is a consequence of the minus sign in the adjoint system (|6]). 

When the bra vector at the initial time is (i^ix-i 0)1 = (1? lj 1> • ■ •)> the generating 
function F(x,t) is expressed in terms of the bra and ket state vectors as follows: 

F(x,t) = (fj( x ,0)mx,t)). (8) 

This fact means that it is sufficient to compare the bra state vector at the initial time, 
(?/>(x,0)|, and the ket state vector at time t, \i/)(x,t)), in order to obtain the generating 
function. Additionally, using the initial conditions for the bra and ket state vectors, and 
recalling the conserved norm (JTJ), we obtain 

(^(x,t)\i>( x ,t)) = (V(x,o)|^(x,o)) = 1. (9) 



3. Cyclic evolution 

If the system develops in a cyclic way, the discussions become easy. For the cyclic case 
with period T, the cyclic state is written as 

|^( X ,T)> = e^V(x,0)>, (10) 

i.e., the ket state vector at time T is proportional to that at the initial time. Note that 
r)(x) is a real number due to the construction of the modified transition matrix H(x, t). 
We here call r)(x) & 'phase'. The reason for this choice of terminology will be explained 
later. 
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(11) 



where we used the conserved norm ([9]). Hence, the phase r](x) is the cumulant generating 
function. As shown in [22H2T] . the cumulant generating function is divided into two 
parts, the dynamical part and a geometric part. The geometric part can be interpreted 
as the Berry phase or Aharonov-Anandan phase for the adiabatic or nonadiabatic 
evolutions, respectively. It has also been clarified that these geometric phases are 
physically-meaningful, and they are deeply related to the pumping phenomena discussed 



We here comment on the reason of calling the quantity f](x) a 'phase'. If one 
considers a Schrodinger equation with Hermitian Hamiltonian, a phase £ is defined 
as exp(iC), where £ is a real number. However, in general cases with non-Hermitian 
Hamiltonian, the phase becomes a complex phase [36] • In this sense, the phase r/(x) is 
considered as a special case of a complex phase. In addition, the dynamical phase and 
geometric phase are adequately defined in a similar way to that in quantum mechanics. 
Hence, the term 'phase' is used. 

4. Geometric phase for noncyclic and nonadiabatic evolution 

For noncyclic cases, we cannot use the geometric phase concept developed in [22H2T] . 
The difficulty basically stems from the fact that the state vector at time t is not 
proportional to that at the initial time. Hence, it is impossible to simply compare 
the 'phase' difference between the two states. 

We here develop a general framework of a nonadiabatic and noncyclic geometric 
phase for counting statistics. After giving some discussions about dynamical and 
remaining phases, we will show the reason why the remaining phase can be interpreted 
as a geometric phase, following closely the discussion given by Samuel and Bhandari [28] . 

Jf.,1, Definition of the remaining phase 

First, we define the dynamical phase as usual: 



The new state vectors, from which the dynamical phase is removed, are defined as 
follows: 



in [22H2Z]. 




(12) 




(13) 



(14) 



As discussed in section 2.2, the norm is conserved and hence we obtain 



^Mx,t)|0(x,t)} = o, (^( x ,tMx,t)> = i. 



(15) 



Note that 5(x, 0) = 0. 
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Using the new state vectors, the generating function can be rewritten as 



Hx,t) = (p(x,0)|exp 



5( X ,t')dt' 



\<f>(x,t)) 



exp 



o 



S( X ,f)dt' 



exp[ln(^(x,O)|0(x,t)}]. (16) 
The total phase fi(x, t), defined as e^ x '^ = -F(x, t), is adequately divided into two parts: 
Kx,t) = [* 5( X ,t')dt' + ]n(tp(x,0)\(f)(x,t)). (17) 
Hence, the remaining phase j(x, t) is defined as 

-Y(x,t) = M<p(x,oMx,t))- (18) 

If the solutions of fl5]) and ([6]) are obtained, it is possible to estimate these two phases. 

Although the two 'phases' were defined, one may wonder whether this division is 
important and meaningful or not. As is easily seen, these phases correspond directly to 
the cumulant generating functions. Hence, one may expect that the remaining phase is 
also meaningful, in a similar way to the cyclic cases in section 3. In the cyclic cases, 
the remaining phase can be interpreted as a geometric phase, and it gives important 
information about the pump current [22H27] . In the next subsection, we will show that 
it is actually possible to interpret the phase 7(x, as a noncyclic geometric phase. 

4-2. Interpretation of the remaining phase as a noncyclic geometric phase 

Let M denote the set of all states |0(x, t)): Af = {|0(x, i)}}- TZ is the space of rays: 
TZ = N I ~, where ~ denotes that elements of A/", which differ only by a phase, are 
regarded as equivalent. In the present paper, the phase is always defined as a real 
number because of the construction of the modified transition matrix. Hence, the 
relation ~ means that an element of M is proportional to another element of M . A 
natural projection map it : M — > 1Z connects each vector to the ray on which it lies. 
The triplet (A/", Tt, 7r) forms a principal fiber bundle over the base space 1Z. Since |0(x, t)) 
satisfies 

^ t \<f>(x,t)) = [H( X ,t) - S(x,t)} |0(x,t)>, (19) 

we obtain the following parallel transport law 

Mx,t)|A|0( x ,t)) = O , (20) 

which defines a natural connection on this fiber bundle. 

Let \4>(x,s)) be a curve in J\f. The tangent vector to this curve is defined 
as \u) = (d/ds)\<f)(x, s)). The corresponding tangent vector in the dual space is 
(v\ = (d/ds)(y?(x, s)\. We here introduce a quantity 

A(x,s) = (<p(x,s)\u( X ,s)). (21) 
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time evolution 







1 
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i geodesic curve r 



K 

\kx,t)) = \<M 



Figure 1. A closed curve for the geometric-phase interpretation. The forward path 
from |^>(xj 0)) to \(j)(x,t)) obeys the time evolution of the generating function, and the 
backward path from |</> ) = \<fi(x,t)) to = \(j)(x,0)) is generated by the geodesic 
curve in 1Z. 



For a cyclic case, the Aharonov-Anandan phase is defined by using A(x, s). In addition, 
introducing gauge transformations, 

\^X,s))^exp[a( X ,s)Mx,s)), (22) 

(<P{X, s )\ -> exp[-a(x, s)](</?(x, s)|, (23) 

the quantity s) transforms as 

da(%, s) 



A( X , s) ->• A(x, s) 



ds 



(24) 



In order to see the geometric characteristics of the phase "f(x, t) defined by ( 1T81) . we 
consider the 'closed' curve in M as shown in figure [TJ That is, the state develops from 
|0(x>O)) to t)) according to the time evolution given by (fl9l) . and the final state 
|0(x, t)) = |0 O ) is connected to the initial state |0(%,O)) = using a geodesic curve. 
The construction of the geodesic curve is as follows. When \<f)(x, s)) is a curve in J\f and 
\u) is its tangent vector, its projection orthogonal to the fiber is 

\Ax, s )) = s )) ~ MX, s)\(p(x, s)) = \u( X , s)) - |0(x, s)){<p(x, s)\u(x, *))■ (25) 

The corresponding projection in the dual space is written as follows: 



(v'(x,s) \ = —(<p(x,s)\ 
as 



{<p( X , s)\A( X , s) = {v( X , s)\ - (<p(x, s)\(v(x, s)\(j)(x, s)). (26) 

Note that |u) and (v\ are not gauge covariant, but \u') and (v'\ transform covariantly 
under the gauge transformation. Because (v'\u') is gauge invariant, we can use it to 
define a metric on 1Z. The metric then determines geodesies in TZ, which are found by 
variation of / (v'\u')dl, where I is an affine parameter determined by the metric. Hence, 
the geodesic equation is given as follows: 

±\u>(x,s)}-A(x,s)\u'( X ,s)) = 0, (27) 
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±(v'( x ,s)\+A(x,s)(v'( X ,s)\=0. (28) 
as 

Using the geodesic curve as defined above, we can show that the phase 7(x,i) is 
written as 

j{ X ,t) = <fA( X ,s)ds. (29) 

In order to show (|2"9~|) . let r(s) be a geodesic curve in 1Z connecting 7r(|0o)) to 
7r(|0i)), and < s < 1. We define the horizontal lift |0 S ) of this curve, starting 
from |0 O ) = l^o) an d we have A(x,s) = for arbitrary s. In addition, a gauge 
transformation |0 S ) = exp[/? s ]|0 s ) is considered, where f3 s is chosen so that f3 = 
and Pi = \n{ip(x, O)|0(x, t)). Since the geodesic equation fl27j) is gauge covariant, \<f) s ) is 
also a geodesic curve, which connects \4> ) and Hence, we have 

l 1 A( X ,s)ds^ NA(x,s) + ^-)ds = | o 1 ^ds = ln(^( X ,O)|0( X ,t)). (30) 

Since the quantity s) is always zero during the time evolution from |0(x, 0)) to 
\4>{Xi we can conclude that the integral in (12"9~|) on the closed line defined in figured] 
adequately gives the phase 7 in f|T8|) . 



5. Concluding remarks 

In the present paper, we gave a general discussion of a noncyclic geometric phase for 
counting statistics, which enables us to extend the adiabatic discussions in j32j|33] to 
nonadiabatic cases. Firstly, we clarified that the counting statistics can be reformulated 
so as to compare a bra state at the initial time and a ket state at an arbitrary 
time. Secondly, removing the dynamical contribution from the cumulant generating 
function, a remaining phase was introduced. The dynamical and remaining phases would 
correspond, respectively, to cumulant generating functions with different meanings. 
Thirdly, it was clarified that the remaining phase is also mathematically meaningful; 
the remaining phase can be interpreted as a 'geometric phase'. In order to interpret the 
remaining phase as the geometric one, a geodesic curve was used, and we obtain a closed 
contour of the integral; the contour of the integral is given by the actual evolution of 
the generating function and a geodesic curve which connects the initial and final states 
of the generating function backwardly. In addition, the remaining phase is clearly gauge 
invariant, and it depends only on the geometric path in the base space; it does not 
depend on its rate of traversal. Hence, the definition of the remaining phase readily 
implies a geometric interpretation. 

We here note that the discussions using the geodesic curves are used only for the 
justification of the geometric property of the phase. In practice, there is no need 
to estimate the geometric phase by finding the geodesic curves; it is enough to use 
(fT8]) . Hence, in order to obtain the dynamical phase and the geometric phase, only the 
solutions of the generating functions are necessary. 
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In section 1, we noted that a mathematically interesting structure may also be 
physically-meaningful. For cyclic cases, as denoted in section 3, the mathematical 
structure is indeed physically-meaningful; the geometric phase gives information about 
the pumping phenomena induced by periodic perturbations [22H27] . It is still unclear 
whether the geometric phase is also important for general nonequilibrium states. In 
order to confirm it, studies for specific models would be useful and needed; these are out 
of the scope of the present paper, and will be clarified in future works. In what follows, 
we will give some discussions to apply the geometric phase to studies of nonequilibrium 
states. 

When adiabatic cases are considered, the dynamical phase gives simply information 
about the number of specific transitions in quasi-static cases. Hence, the noncyclic 
geometric phase will correspond to an additional contribution from the time- dependent 
transition matrix. Actually, in j32j|33], such effects in adiabatic cases are discussed for 
a simple two-state model. This means that even in a quasi-static evolution, the total 
number of transitions is not given by a quasi-static contribution. It would sometimes 
be beneficial to discuss these two contributions separately, and the noncyclic geometric 
phase discussions divide them naturally even in nonadiabatic cases. 

Next, let us consider transitions between two nonequilibrium steady states. In such 
cases, it has been pointed out that a concept of 'excess' quantities would be important to 
construct useful identities or relations; excess heat or entropy production have already 
been discussed in a few contexts [13 | [T ^l371 - H2] . Here, the 'excess' quantity means an 
additional quantity caused by time-dependent perturbations or external fields, and it 
gives no contribution if the state is a steady one. For example, an excess heat has been 
explicitly introduced and calculated in a Langevin system, and the concept is deeply 
related to the Hatano-Sasa relation [13] , which can be considered as an extension of the 
Jarzynski equality to nonequilibrium cases. The excess heat is estimated by subtracting 
a steady heat, the so-called housekeeping heat, from the total heat. The housekeeping 
heat is needed to maintain the steady state out of equilibrium. One of the remarkable 
characteristics of the excess heat is that it is a finite quantity. In contrast, the total and 
the housekeeping heat during time t diverge when t — > oo. It may be plausible to use 
such finite quantities in order to discuss nonequilibrium steady states j37j. An important 
point here is that a total quantity is divided into two parts, i.e., the housekeeping one 
and the excess one. In the noncyclic geometric phase discussions, we observed that the 
total phase is naturally divided into two parts, i.e., the dynamical one and the geometric 
one, as described above. In this sense, one may expect that this natural division of the 
total phase have some physical meanings, as in the case of the excess quantity. Although 
further studies will be needed to clarify the physical meanings of the noncyclic geometric 
phase, we hope that such discussions based on the noncyclic geometric phase would give 
deep insights for nonequilibrium states. 
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Appendix A. Generating function for counting statistics 

Although counting statistics has been discussed in various contexts (for example in [T9]). 
we here give a rederivation of the generating function used in counting statistics for the 
reader's convenience. 

In the framework of counting statistics, the quantity of interest is the number of 
target transitions. The target transitions can be chosen arbitrarily. For example, if 
we consider a difference between the number of transitions from state i to state j (we 
denote this transition as i — > j) and that of transition j — > i we can calculate a net 
'current' of i — > j. The information about such 'current' plays an important role in 
discussing nonequilibrium states. 

Let {K nm (t)} be a time-dependent transition matrix. Although one can calculate 
the statistics for multiple target transitions in general, we here derive the generating 
function for counting the number of events of a specific target transition i A — > Ja, for 
simplicity. The generalization to multiple transitions is straightforward [19J. 

Let us denote the probability, with which the system starts from state m and 
finishes in state n with Aa being the number of target transitions i A — > 3a during 
time t, as P nm {^A\t)- In order to calculate the probability P n m{NA\t), we here define 
a probability G' kl (t) with which the system evolves from state / to state k, provided no 
3a transitions occur during time t. By using the probability G' kl (t), the probability 
Pnm (Aa|£) is calculated as 

Pnm(N A \t) = G' njA (t) * K jAiA (t)G' iAjA (t) * • • • * K jAiA (t)G[ AjA (t) *K jAiA (t)G' iAm (t), (A.l) 

jV A -l 

where gi(t) * g 2 (t) = J ' g\(t — t')g2(t')dt' denotes the convolution. This formulation 
means that an occurrence of the target transition i^ — > jx is sandwiched in between 
situations with no occurrence of the target transition, and it is repeated Aa times. 
Next, we construct the generating function f nm {x^) °f t ne probability P nm (AA|t): 

oo 

fnm(x,t)= E e^P nm (N A \t). (A.2) 

N A =0 

That is, the generating function f nm (x, t) gives the statistics of the number of transition 
^a - > 3a during time t under the condition that the system starts from state m and ends 
in state n. The generating function f n m{Xi^) satisfies the following integral equation 

f nm (x, t) = G' nrn (t) + f G' njA (t - t')e*K jAiA (t')f iAm ( X , 0< (A-3) 
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and obeys the following time-evolution equation 
d 



d fnm(x,t) =J2K ni (t)GUt) - S n>jA K jAiA (t)G[ Am (t) + eXG> njA (0)K jAiA (t)f iAm (t) 



+ f Q (i^At-t')) e*K jAiA (ty~ iAm (t')dt' 
= Y, K m(t)fim{x,t) - <W* - e*)K jAiA (t)f lAm (x,t), (A.4) 

i 

where f nm (x, 0) = 5 n m . In order to show (IA.4[) . we used the following two facts: Firstly, 
the probability of no target transitions, G' nm (t), obeys 

^G' nm (t) = K ni (t)G' im (t) - S n , jA K jAiA (t)G' iAm (t), (A.5) 

where G' nm (0) = 5 nym . Secondly, the derivative of the convolution is given by 

^ J^g^t - t')g 2 (t')dt' = 9l (0)g 2 (t) + jf* (^ 9l (t - A g 2 (t')dt>. (A.6) 

Using the generating function /„ m (x,t), we construct the restricted generating 
function f n (x^) i n © an d © as follows: 

fn(X,t)=J2fnm(X,t)Pm(0), (A.7) 
m 

where p m (0) is a probability distribution at initial time t = 0. Hence the generating 
function is given by ©. In addition, from ( 1A.4|) and ( IA.7j) . the restricted 

generating function f n {x^) satisfies 

f t fn(x,t) =J2 K m(t)fi(x,t) - S n>jA (l - e*)K jAiA (t)f iA ( X ,t), (A.8) 
and these equations should be solved with initial conditions f n (x^) — 

I2 m fnm(X,Q)Pm(0) =Pn(0). 
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